Fundamental Structural Constraint of Random Scale-Free Networks 
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We study the structural constraint of random scale-free networks that determines possible com- 
binations of the degree exponent 7 and the upper cutoff k c in the thermodynamic limit. We employ 
the framework of graphicality transitions proposed by [Del Genio and co-workers, Phys. Rev. Lett. 
107, 178701 (2011)], while making it more rigorous and applicable to general values of k c . Using 
the graphicality criterion, we show that the upper cutoff must be lower than k c ~ N 1 ^ for 7 < 2, 
whereas any upper cutoff is allowed for 7 > 2. This result is also numerically verified by both the 
random and deterministic sampling of degree sequences. 
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Complex networks [l[ are found in diverse natural and 
artificial systems, which consist of heterogeneous ele- 
ments (nodes) coupled by connections (links) markedly 
different from those of ordinary lattices. In particular, 
many systems 2- 6] can be interpreted as scale- free net- 
works, in which the fraction of nodes with degree k (i.e., k 
links) obeys the power-law distribution P(k) ~ /c~ 7 over 
a broad range of values bounded by k m < k < k c where 
7 is called the degree exponent, k m is the lower cutoff, 
and k c is the upper cutoff. There have been interests 
in topological and dynamical properties induced by the 
degree distribution, which have been examined through 
various studies on random scale-free networks 

Random scale-free networks refer to an ensemble of 
networks constrained only by the parameters 7, k m , and 
k c . In general, k m is set as a constant, while k c is assumed 
to increase with the number of nodes N as k c ~ N a 
with Q < a < 1. Besides, self- loops or multiple links 
between a pair of nodes are often disallowed. Under the 
circumstances, the degree exponent 7 and the cutoff ex- 
ponent a determine various properties of networks in the 
thermodynamic limit, N — > 00. It is known that 7 con- 
tributes to the resilience against node failures [7|, the 
epidemic threshold Q , the consensus time of opinion dy- 
namics etc. Meanwhile, a affects the expected value 
of the generated maximum degree [l(J, degree correla- 



tions 



ll| , finite-size scaling at criticality [12], ll3| > etc. 



The studies of random scale-free networks character- 
ized by 7 and a must be based on the knowledge that 
such networks actually exist in the thermodynamic limit. 
Hence, it is necessary to understand the constraint on the 
possible values of 7 and a. This problem is exactly equiv- 
alent to the issue of the graphicality of random scale-free 
networks. A degree sequence {k\, k%, . . . , k^} is said to 
be graphical if it can be realized as a network without 
self-loops or multiple links. As an indicator of the ex- 
istence of graphical sequences, the graphicality fraction 



g [lj] is defined as the fraction of graphical sequences 
among the sequences with an even degree sum generated 
by P(k). Note that the degree sequences with an odd 
degree sum are left out, since such sequences are trivially 
nongraphical. The constraint on the possible values of 
7 and a can be obtained from the behavior of g due to 
the fact that the random scale-free networks with given 
7 and a exist in the thermodynamic limit if and only if 
g is nonzero as N — > 00. Using the graphicality criterion 
given by the Erdos-Gallai (EG) theorem [l5| , Del Genio 
and co-workers 14| have recently studied the behavior of 
g as a function of 7 only for the special case of k c = N— 1 
and k rn = 1. They found 



9 = 



if < 7 < 2 

1 otherwise, 



where the abrupt changes of g at 7 = and 7 = 2 were 
termed graphicality transitions [LJ]. This result implies 
that there exist only sparse random scale-free networks 
with a finite average degree (7 > 2) and left-skewed net- 
works with an abundance of hubs (7 < 0) in the thermo- 
dynamic limit when the range of degree is kept maximal. 

In this Letter, we generalize their study to arbitrary 
choices of degree cutoffs, so that we can provide the com- 
plete picture of the constraint on the possible values of 7 
and a. Starting from the EG theorem, we present a rig- 
orous derivation of g as a function of both 7 and a, which 
is then verified and supplemented by numerical results. 

The EG theorem [l5[ states that a degree sequence 
sorted in the decreasing order k\ > ki > . . . > k^ is 
graphical if it has an even sum and satisfies the EG in- 
equalities given by 



JV 



ki < n(n — 1) + min [n, fc,] 



(1) 
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for any integer n in the range 1 < n < N — 1. 
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TABLE I: Scalings with N of the nth largest degree k n for 
arbitrary values of a, j3, and 7. 
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To determine whether those inequalities are satisfied 
in the thermodynamic limit, we derive the network-size 
scalings of the left-hand side (lhs) and the right-hand side 
(rhs) of each inequality. As the first step, we calculate 
the scaling of k n from the cumulative mass function of k n , 
which is denoted by Tm ■ The maximum degree (n = 1) 
satisfies 
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r^(fc) = n Prob [ki<k] = [C(k)\ 



N 



(2) 



where C is the cumulative mass function of degree. Since 
satisfies the recursive relation 
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(k) = Prob [k n -\ > k and k n < k] 
N-n+l 



(3) 



we can obtain its exact form as 

n-l 



r^(fc) = 



(4) 



Suppose n = i/jV* 3 , where v > and < (3 < 1. If /3 > 0, 
we can use the following approximation for large N: 
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(5) 



Using Eq. Q for /3 = and Eq. © for (3 > 0, we can 



find the rang eof k in which (A:) increases from to 1 



TABLE II: Scalings with N of each side of the nth EG in- 
equality for arbitrary values of a, f3, and 7. 
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in the limit N — >• 00. Since the typical values of fc n must 
fall within this range of k, we can obtain the network-size 
scalings of k n as listed in Table Q] 

While both sides of the EG inequalities are sums over 
k n , we can approximate those sums as integrals, since it 
does not affect the the leading iV-dependent term that 
determines the scaling relation. It is straightforward to 
approximate the lhs, while the rhs needs a careful refor- 
mulation. The second term of the rhs satisfies 
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min [n, k t ] = NO (n - k m ) ^ kP ( k ) 
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+ N0 {k n+1 - n - 1) ( 6 ) 

max[n+l,/j m ] 

where 9 denotes the Heaviside step function defined by 
6(x) = 1 if x > 0, and 6(x] = otherwise. From now on, 
each summation can be converted to an integral over the 
same range. Calculating all the integrals, we can single 
out the leading iV-dependent terms of each side, as listed 
in Table HU The scalings of those terms are completely 
determined by the three exponents a, /3, and 7, while the 
lower cutoff k m turns out to be irrelevant. 

We can now determine whether the EG inequalities are 
satisfied through the comparison of scaling exponents in 
both sides. By the EG theorem, g = 1 if the inequalities 
are satisfied for all possible values of /3, and g = if 
there exist the values of f3 at which some inequalities are 
violated. Hence, the asymptotic behavior of g is obtained 
as follows: 



ifl/a< 7 <2 

1 if 7 > 2 or a < min [I/7, 1] 



(7) 



We note that the behavior of g for the special case of 
a = 1 and 7 < 1 cannot be determined by our scaling ar- 
gument, since both sides of the EG inequalities satisfy the 
same network-size scalings. To address this problem an- 
alytically, it is necessary that we consider the coefficients 
of the leading TV-dependent terms, which is beyond the 
scope of this Letter. Instead, we settle for its numerical 
resolution at the end of this Letter. 

For the other cases, we can analytically determine the 
locations of graphicality transitions from Eq. ([7]). There 
exist two transition points for each value of a in the range 
1/2 < a < 1, namely the upper transition point 7^ = 2 
and the lower transition point 7^ = 1/a. On the other 
hand, no transition occurs for < a < 1/2 where g = 1 
always holds. 

All the predictions on the asymptotic behavior of g 
can be numerically checked by the evaluation of the EG 
inequalities. First of all, we can verify all the scalings 
listed in Table HI some of which are shown in Fig. Q] This 
indirectly supports our predictions on the behavior of g, 
as all the predictions were deduced from those scalings. 
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To obtain direct support for our predictions, we need 
to measure the 7 dependence of g from the random sam- 
ples of degree sequences, as illustrated in Fig. [2] Due 
to sample-to-sample fluctuations at finite system size, g 
changes continuously between and f over a finite range 
of 7, which becomes narrower as N increases. In fact, g 
may not even reach if N is too small, as exemplified 
by the curves for a — 0.825 in Fig. [2ja). The curves in 
Fig. EJb) show that the minimum of g gradually reaches 
zero as N increases. Keeping those observations in mind, 
for the sake of convenience, we regard the range of 7 in 
which g falls below 0.99 as the effectively nongraphical 
region. Then, the boundary of this region can be chosen 
as the effective transition points at finite N, which are 
again marked as 7£ and 7u in Fig. [2fa). 

We also consider deterministically generated degree 
sequences defined by I — C(k n ) — v, obtained from 
Eq. (0, which ensures that the sequences exactly fol- 
low the network-size scalings of k n for n — vN . Those 
sequences filter out the sample-to-sample fluctuations, 
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FIG. 2: (Color online) 7 dependence of the graphicality frac- 
tion g obtained from 10 3 random degree sequences, (a) a 
increases by steps of 0.05 with N = 2 4 x 10 3 . 7£ and 7u 
indicate the two transition points at a = 0.825, whose N 
dependence is shown in (b) as N increases by factors of 2. 
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FIG. 1: (Color online) N dependence of the nth EG inequality 
obtained from 10 3 random degree sequences (symbols) for 7 = 
1.5 and a — 0.9. It is in good agreement with the scalings 
(lines) predicted in Table [II] which is also confirmed by the 
insets: The ratio between the successive slopes of each symbol 
and the slope of the corresponding line stays close to one. 



making it very straightforward to locate the effective 
transition points [see the inset Fig. |3[a)]. They also 
greatly improve the efficiency of calculation, allowing us 
to check our predictions at larger N. We observe that the 
transition points estimated by both randomly and deter- 
ministically generated degree sequences approach each 
other as N — > 00 (for example, see Fig.0|). Therefore, we 
can use either of those two different samplings to numer- 
ically check our predictions. 

In Fig. [3l we present graphicality diagrams obtained 
at two different values of the cutoff coefficient c, where 
the transition lines at finite network sizes are estimated 
by the deterministic samplings of degree sequences, and 
also compared with the transition lines in the thermo- 
dynamic limit predicted by the scaling argument. The 
numerically estimated transition lines tend to approach 
analytically predicted ones, which are independent of c. 
Combining this observation with the verification of the 
scaling relations listed in Table [TTl we can safely conclude 
that the numerical results are slowly converging to our 
predictions as N increases. 

Moreover, Fig. Ogives us some clues as to the location 
of the graphicality transitions for a — 1 and 7 < 1, which 
could not be determined by the scaling argument as pre- 
viously explained. The estimated transition lines suggest 
that the location of the transition is dependent on c : 7* 
approaches 7 = at c = 1, as previously reported [l4|, 
but it converges to some different limiting value if c =/= 1 . 
The effect of c on the graphicality at a = 1 is more closely 
examined in Fig.[4ja), which suggests that 7^ varies con- 
tinuously between and 1 with c, while 7^ converges to 
2, regardless of c. This nicely contrasts with Fig. HJb), 
which confirms our prediction that the transition points 
are independent of c for a < 1. 
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While we have given an almost complete picture of 
the graphicality issue of scale-free networks, the nature 
of graphicality transitions requires further studies. Note 
that at a transition point, the comparison of scalings fails 
to determine whether the Erdos-Gallai inequality holds 
in the asymptotic limit, just like the case of a = 1. In 
such cases, the coefficients of the leading-order terms as 
well as the second- and higher-order terms must be con- 
sidered to determine the value of g. Thus, we cannot 
claim yet that graphicality transitions are truly discon- 
tinuous as previously claimed [l4j |. since they might be 
sharp but continuous transitions resembling the contin- 
uous change of 7£ with c at a = 1. The claim should 
be either proven or disproven by a more complete under- 
standing of the behavior of g at the transition points. 

In conclusion, we have found that in the thermody- 
namic limit random scale- free networks without self-loops 
or multiple links are either sparse (7 > 2) with arbitrary 
values of degree cutoffs, or dense (0 < 7 < 2) with the 
upper cutoff k c ~ N a satisfying a < I/7, supplement- 
ing the statement that "all (random) scale free networks 
(with maximal range of degree) are sparse." 14J This 
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FIG. 4: (Color online) N dependence of the transition points 
as c increases by steps of 0.2 for (a) k c = c(N — 1) and (b) 
k c = ciV ' 95 . Symbols indicate the transition points obtained 
from up to 10 3 random degree sequences, while those obtained 
from deterministic sequences are connected with lines. 



also a gre es with the upper cutoff found by Seyed-allaei 
et al. [16], which is required for scale-free networks with 
7 < 2 generated using a node- fitness mechanism 17f ■ We 
also numerically found that the cutoff coefficient c affects 
the realizability of degree sequences for the special case 
of the linear cutoff a = 1, which has been overlooked. 
Our results impose a limit on the values of 7 and a for 
which the properties of random scale-free networks nu- 
merically obtained in finite systems can be extrapolated 
to the thermodynamic limit. 
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FIG. 3: (Color online) Graphicality diagram in the (7, a) 
plane at (a) c = 1 and (b) c = 0.6. The nongraphical region is 
in gray, while the other lines are the transition lines estimated 
using deterministic degree sequences, whose behavior of g is 
as illustrated for a = 0.95 in the inset of (a). 
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